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Proced6 de localisation spatio-temporelle, de radio-mobiles en milieu urbain 


(57) 


Le proc6d6 consiste : 





k calculer (1) la r^ponse impulsionnelle du canal de 
propagation entre te radio-mobile et la station 6met- 
trice/r6ceptrice, 

k estimer (2) a partir de la r6ponse impulsionnelle 
les couples (angles, retard) du trajet de propagation 
direct et des multitrajets r^sultats de reflexions sur 
des obstacles situ£s dans le champ de propagation 
des ondes entre la station de base et le radio- 
mobile, 

& calculer (2) pour chaque trajet les attenuations otj, 
et & calculer (3) la distance d sSparant la station 
6mettrice/r6ceptrice du radio-mobile en consid6- 
rant seulement les retards et les attenuations du 
trajet direct et d'un multitrajet quelconque. 

Applications : Localisation et suivi de radio-mobiles. 
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Description 

[0001] La pr6sente invention concerne un proc6d6 de localisation de radio-mobiles en milieu urbain. 
[0002] Elle s'applique notamment k la localisation et au suh/i de radio-mobiles dans un r6seau GSM pour anticiper 
5 leur passage d'une cellule a la suivante et ameiiorer la reception des stations de base par la connaissance de la posi- 
tion des radio-mobiles 6metteurs. 

[0003] Dans des configurations simples de propagation, la localisation de radio-mobiles est traditionnellement obte- 
nue a I'aide de radio-goniometres disposes en des lieux determines d'un terrrtoire. La position d*un radio-mobile peut 
alors §tre d6termin6e par triangulation en reportant sur une carte geographique les lieux oil sont situes les radiogonio- 

10 metres ainsi que les directions d'arriv6e sur ceux ci des ondes radio6lectriques 6mises par le radio-mobile. 

[0004] Cependant, cette methode qui repose sur ['exploitation des seuls chemins de propagation directs apparaTt drf- 
f icilement applicable dans le cas ou les ondes empruntent plusieurs trajets et est egalement lourde de mise en oeuvre 
car elle n6cessite plusieurs stations de reception. Ce type de propagation est notamment celui qui est rencontre en 
milieu urbain dans les reseaux de communication radio-mobiles cellulaires, ou le spectre des communications est situe 

75 dans une gamme d'ultra hautes frequences s'etendant typiquement de 30 MHz a 3 GHz pour le reseau GSM. Pour 
cette gamme de frequences la propagation des ondes radioeiectriques entre un emetteur et un recepteur peut 
empreinter des trajets multiples dus aux reflexions et aux diffractions sur des obstacles interposes sur le chemin de pro- 
pagation direct entre remetteur et le r6cepteur. La gamme des frequences utilises est aussi la cause de trajets drffus, 
aussi bien temporellemerrt que spatialement, qui degradent les performances des liaisons. 

20 [0005] Le but de I'invention est de pallier les inconvenients precites. 

[0006] A cet effet, I'invention a pour objet un proc£d6 de localisation spatio-temporelle d'un radio-mobile relativement 
a une station de base emettrice/receptrice caracterise en ce qu'il consiste : 

a calculer la r6ponse impulsionnelle du canal de propagation entre le radio-mobile et la station emettrice/r6cep- 
25 trice. 

- . k estimer a partir de la reponse impulsionnelle les couples (angles, retard) du trajet de propagation direct et des 
multitrajets resultats de reflexions sur des obstacles situ6s dans le champ de propagation des ondes entre la sta- 
tion de base et le radio-mobile, 

k calculer I'attenuation q de I'onde propagee par chaque trajet, 
30 - et k calculer la distance d separant la station emettrice/receptrice du radio-mobile en considerant seulement les 
retards et les attenuations du trajet direct et d'un trajet non direct quelconque. 

[0007] D'autres caracteristiques et avantages de I'invention apparaTtront k I'aide de la description qui suit faite en 
regard des dessins annexes qui represented : 

35 

la figure 1 le precede selon I'invention sous la forme d'un organigramme, 
la figure 2, un modeie de reception monocapteur, 

la figure 3, une illustration d'une reception multicapteur par une antenne circulaire, 

les figures 4 et 5, une illustration d'une transmission k deux ou trois trajets en milieu urbain. 

40 

[0008] Le procedS selon (Invention se d6roule de la facon representee k la figure 1 suivant trois etapes principales. 
Dans une premiere etape, representee en 1 k la figure 1 , il consiste k calculer la reponse impulsionnelle h(t) du canal 
de propagation pour estimer a I'etape 2 les angles e m et les retards x m des differents trajets ainsi que les attenuations 
Op des ondes propagees sur chacun des trajets et enf in a calculer k I'etape 3, k partir des retards i m calcuies et de la 
4$ reponse impulsionnelle h(t) du canal, la distance d separant un radio-mobile de la station d'emission/reception. 

[0009] Le calcul de la reponse impulsionnelle du canal qui est effectu6 & I'etape 1 a lieu sur le signal obtenu k la sortie 
d'une antenne multi-capteurs. 

[001 0] En presence de plusieurs trajets M d'emission bande etroite de radio-mobiles, le signal obtenu a la sortie d'un 
capteur n de I'antenne multi-capteurs est 6gal k la somme des M signaux s m (t) 6mis dans les directions 6 m retard6s de 
so x m et multiplies par la reponse aj(e m ) du capteur i dans les directions e m a laquelle s'ajoute un bruit bj(t). Dans I'hypo- 
these d'une transmission numerique a bande etroite de frequence, une representation de la sortie vectorielle de 
I'antenne compos6e de N capteurs s'obtient en concatenant la sortie X;(t) de chaque capteur, chaque sortie Xj(t) v6ri- 
fiant la relation : 

55 M 

*iW= Esjtj.a^ej+b,© (1) 

m=1 



2 
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[001 1 ] En introduisant un vecteur directionnel a(9m) tel que 

a(9 m ) = [a 1 (e m ) a N (e m )]' (2) 

et en concatenant la sortie Xj(t) de chaque capteur, on obtient un vecteur x(t) tel que : 

M 

x(t) = [x,©, Xi(t). x N (t)] 1 = £ a(8Js m (t)+b(t) (3) 

m=1 

[001 2] En rassemblant les vecteurs directionnels £(9,^ dans une matrice 

A=[a(9 1 ) l a(9 2 ) 1(6 n)] 

!a relation pr6c6dente s'6crit encore : 

x(t) = As(t)+b(t) (4) 



avec s(t) = [s 1 (t).s 2 (t)- * -s M (t)]\ 
ou $1 est le signal du 1 er trajet, 
s 2 est le signal du 2 i6me trajet, 

s M est le signal du M'* mG trajet. 



[0013] Suivant le modfeie de propagation qui est represents schematiquement a la figure 2, le signal x(t) qui est 
obtenu a la sortie d'un capteur de reception est un signal de la forme : 



M 

x(t)= £a r 
m=1 



e'^.aOJ.sCt-xJ + bW 



(5) 



ou designe ('attenuation du m l6me trajet 
T m est le retard de groupe du m i6me trajet 
M est le nombre de trajets 
et c|> m est le dephasage du m i6me trajet. 

soit encore sous la forme d'une convolution du signal 6mis s(t) avec la reponse impulsionnelle h(t) du canal : 

soitx n (t) = h n (t)*s(t) + b n (t) (6) 



ou b n (t) reprSsente I'enveloppe complexe du bruit additif. 
[001 4] Dans le cas d'une reception sur I'antenne compos£e de multi-capteurs, comme cela est represents a la figure 
3, la reponse impulsionnelle du canal est modifiee par la reponse de I'antenne suivant les directions d*arriv6e. Cette 
reponse est de la forme : 



h n (t)= £a n (e fTV y m )a m e^6(t-T m ) 



(7) 
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en supposant le canal id£al, c'est-a-dire de bande infinie. 
[0015] En mode 6chanti!lonn6, la nSponse du canal de bande B prend la forme : 

M 

5 h„M»F e £a n (eja m e'* m sinc((k-F e T m )B) (8) 

m=»1 



avec k e Z, I'espace des nombres entiers. 

10 

sin(n(k-FeT m )B) 



etsinc((k-F e , m )B) = n(k-F e T m )B « 



rs F e est la frequence d'Schantillonnage 

B est la bande en frequence r£duite du signal. 

[0016] Cette r6ponse peut £tre estimde par un filtre causal, a niponse impulsionnelle finie de longueur L et de coef- 
ficients tels que : 

20 

M it 

h n [k]= £a n (eja m e* m sinc((k-F e x m )B) (10) 

25 

pourke {0 ... L} 

et h n [k] = 0 pour k < 0 ou k > L 



30 



Le vecteur h n correspondant s*6crit h n = [h n (0), • • • , h n (L)] 1 



(11) 



[0017] Suivant un premier mode de realisation de I'invention, ("estimation des angles e m et des retards x m de l'6tape 
2, est ex6cut6e de fa$on disjointe, en utilisant la relation matricielle suivante : 



35 



H = [h 1 .h 2 ,...,h N ] = M.n.A t (12) 



avec : 



40 



45 



M = 



sinc(-t-|F e B) sinc(-t 2 F e B) 
sinc^l-T^ejB) sinc((1- t 2 F 6 )B) 



sinc(-T M F e B) 
sinc((1-T M F e )B) 



sincaL-x^eJB) sincaL-^Fe) B ) sinc((L-T M F e )B) 



(13) 



50 



55 



4 



n = 
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0 

0 a 2 



: ••. 0 
0 ••• 0 a N 

A=[a 1 (e 1 )...a N (6 N )] 



10 



(14) 



(15) 



[0018] L'estimation des paramfctres temporels x m s'obtient k partir de la relation (12) en estimant un projecteur 



15 



HH 1 



sur I'espace temporel HH* et en appliquant soit la m6thode connue sous l'abr6viation anglo-saxonne MUSIC, 
de < <Multi Signal Classification ))en cherchant les maxima de la quantity 



20 



25 



30 



35 



1 



n?(T) t P HHt m(T) 



(16) 



oil 



m(T) = [sinc(-xF e B) sinc((L - tF c )B)] 1 

est une colonne de la matrice M(13), soit la methode deterministe connue sous le nom de << maximum de vraisem- 
blance ) ) en utilisant un critdre de la forme 

J(tl T M ) = t r ((l L+1 -Wi(MtWi)- 1 M^HH t ) . (17) 



40 



45 



[001 9] Une description de la methode MUSIC peut §tre trouvee dans Particle de M. Ralph O. Schmitt ayant pour litre 
< { Multiple source DF signal processing : An experimental system >> publiee dans IEEE transaction on antennas and 
propagation AP34 N° 3 March 1986. 

[0020] La determination des M couples (e m , i m ) permet de localiser les radio-mobiles. Dans le cas ou, comme cela 
est repr6sent6 k la figure 4, seulement deux trajets, un trajet direct et un trajet comportant une reflexion sont seiection- 
n§s, la determination des couples a lieu en considerant quatre parametres 9 1t t-,, 6 2 , Q u ' se ram^nent k la connais- 
sance de e lf 8 2 et des retards , t 2 . En appelant d la distance du trajet direct qui separe remetteur radio-mobile 4 d'un 
recepteur 5 et , d 2 les distances du trajet indirect separant respectivement le recepteur 5 d'un obstacle 6 sur lequel 
I'onde indirecte est r§fl6chie et remetteur 4 de I'obstacle 6. la distance d est fournie par la relation : 



so 



55 



d= — r 



(t 2 -t,)c 



(18) 



^1-1 

a 2 \i 



ou a-| et cc2 designent les attenuations du premier et du deuxieme trajet et \i la perte due k la reflexion du second 
trajet. Les attenuations a-, et sont obtenues k partir des couples (0 1( t-j ; e 2 , t 2 ) et de la relation matricielle 



5 



r 
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n = M # H A #t = 




(19) 



5 



oCi a,- correspond & 1'attenuation du r mo trajet. 
[0021 ] Comme la perte de puissance dans une propagation libra se fait en 



w 



15 



ou r est la distance parcourue, et a 2 sont Ii6s par la relation : 



20 



«i _ (d 1 4d 2 ) g 

«2 d 2 



(20) 



avecd 



(x 2 -T 1 ) 2 .c 2 +2d(T 2 -x 1 )c 
2[(t 2 * x 1 )c+d-dcos{6 2 -6 ! )] 



(21) 



25 



etd 2 = V d2+d i -2dd 1 cos(e 2 - e n ) 



(22) 



[0022] S'il y a plus de deux trajets, la figure 5 montre que 1'information apportee par tout trajet suppl6mentaire se 
30 resume k deux parametres, et les nouvelles inconnues sont 6galement au nombre de deux. En consequence, une 
detection de plus de deux trajets n'a pas d'interet dans le sens ou elle ne supprime pas les inconnues. II suff it dans ces 
conditions de limiter les traitements en considSrant seulement le trajet direct et un trajet autre quelconque. Connaissant 
les angles d'arriv£e des diff erents trajets. il est possible, grace & une carte topographique d£taiII6e en trois dimensions 
de trouver les points de reflexion les plus probables et ainsi de remonter & la source demission. 
35 [0023] D'autres modes d'estimation des parametres sont egalement possibles. 

[0024] L'estimation disjointe des parametres decrite pr6c6demment a I'avantage de reduire I'espace de recherche. 
Cependant, dans cette methode , les angles et les retards ne sont pas associes si bien qu'elle souffre d'un manque de 
resolution pour estimer les retards proches. 

[0025] D'autres modes de mise en oeuvre de ['invention permettent d'associer les parametres en estimant la matrice 
40 des attenuations a avec les differentes combinaisons (e m> x m ) possibles en cherchant les couples rendant la matrice 
Cl la plus diagonale possible. 

[0026] Selon un deuxfeme mode de mise en oeuvre de I'invention, le procede opere par estimations successives en 
estimant successivement les angles e et les retards t tout en les associant. Pour cela, deux methodes peuvent §tre uti- 
lisees. La premiere consiste & estimer les parametres temporels -t m connaissant les angles 6 m . 
45 [0027] En estimant des angles 6 m par la methode MUSIC ou la methode du maximum de vraisemblance l'estimation 
des retards t m est obtenue en estimant la matrice M par la relation matricielle 



[0028] Chaque colonne de la matrice M forme un projecteur permettant de trouver le retard correspondant x m en 
appliquant I'algorithme MUSIC pour une seule source d'6mission. 
55 [0029] Le critere utilise est de la forme : 





50 
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P(t) = 



Btft^Pi, m(t) 
m^(-t) n](x) 



(26) 



10 



ou P - - gK^i)03 T (Ti) 
^(TiJrrKtj) 
avec m(x) tel que M = [rn(xi)... rn(tM)] 
sinc(-x F e B) 



et m(T) = 



15 



sinc((L-xF e )B) 



[0030] L' estimation des retards peut 6galement etre obtenue par correlation. En considerant que le signal recu sur 
un capteur de M sources est de la forme : 



20 



25 



30 



35 



40 



M 



x(t)= £a(e m )s m (t) + b(t) 

m=1 



et connaissant les angles 6 m , il est possible de considerer un signal : 

s[k] = A # x[k] 

qui est forme de M versions du signal 6mis retardees tel que : 



s(k-x 2 F e ) 
s(k-T M F e ) 



(26) 



(25) 



[0031] (-'autocorrelation d'une des composantes de s [k] permet d'obtenir une reference temporelle et I'intercorrela- 
tion de cette composante avec les autres composantes donne une estimation des retards relatifs. 
45 [0032] De fagon similaire a la methode precedente, il est aussi possible d'estimer les parametres spatiaux e m con- 
naissant les retards t m . Ceux-ci sont obtenus en calculant I'equation matricielle : 



so 



A 1 =(M f M) 1 MtH 



(27) 



[0033] Chacune des melhodes utilisees dans ce deuxieme mode de realisation du procede seion invention permet 
apres avoir estime les angles ou les retards de trouver les autres inconnues par M recherches en une dimension. Elles 
55 ont pour avantage de repousser certaines ambiguTtes, mais la seconde recherche apparaTt tres d6pendante de fa qua- 
lity de la premiere estimation. Une alternative serait d'effectuer, comme cela a 6te d§crit par M. J. GROUFFAUD dans 
son m6moire de th6se de Doctorat de I'Ecole Normale Sup6rieure de Cachan de juin 1997 intrtul6e <( Identification 
spatio-temporelle de canaux de propagation a trajets multiples ) ) , des recherches successives en alternant les m6tho- 



7 
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des pour aff iner la precision de i'estimation des parametres. 

[0034] Enfin, suivant un troisieme mode de realisation du proc&ie selon Hnvention, I'estimation des parametres 9 m 
et T m a lieu de fa^on conjointe en appliquant des methodes de determination des reponses impulsionnelles fondees sur 
une decomposition en sous-espaces, un sous-espace signal de dimension M egale le nombre de sources et un sous- 
espace bruit complementaire du sous-espace signal de dimension N(L + 1)-M, ou en appliquant la methode du maxi- 
mum de vraisemblance. 

[0035] Pour estimer conjointement tes angles et les retards par la methode des sous-espaces, un project eur est 
cree a partir de la matrice d'autocorrelation des f iltres 



10 



15 



20 



[0036] L'estimation 6 du vecteur h = [h ^O], h ,[1], .... h^L] h N [0], h N [1] h u [L]] x a lieu en recherchant les 

maxima de la quantity 



p (h) = — ! (28) 

(h) gt(0.x).P Rh .g(e,t) 



ou g(e,T) = a(e)®m(*0 et 

25 M 



b= £ q m .g(e m ,T r 

m=1 



avec 



[0037] Dans la recherche des parametres par la methode du maximum de vraisemblance les parametres cherches 
35 sont les minimum du critere, 

J (Bi. T! 6 m . x m ) = tr [l^ L+i) - G(G+G)" 1 Gt R h J 

ou G = [g(9 1l t 1 )-g(e m> T m )] 



45 [0038] Cette recherche a lieu dans Tespace de dimension M x 2 des parametres {Q^ , tj . * • • , e M , x^. 

[0039] L'optimisation des differents modes de realisation du proc6de selon I'invention peut avoir lieu en appliquant 
une methode connue d'optimisation du type Gauss-Newton. Apres avoir trouve un point d'initialisation par les methodes 
decrites precedentes, une optimisation des angles 6 m peut avoir lieu en deroulant I'algorithme suivant : 

so Initialisation : 6 0 = e injt 
Iteration k : 

Calcul du gradient v k .-j et du hessien au point e k .-j 
Determiner 



e k =e k _ 1 -H-: 1 v k _<, 
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Test de la convergence 

si |V k _ 1 | <e STOP anon k = k + 1 

5 

[0040] Le m£me type d'optimisation peut etre utilise pour optimiser la recherche des retards t m . 
[0041] Un autre mode d'optimisation peut §tre mis en oeuvre par une version approch§e du maximum de vraisem- 
blable suivant par exemple la m 6th ode d6crite dans IEEE transaction on signal processing, Vol. 39, N1 1, Novembre 
1 991 des auteurs M. Viberg, B. Ottersten, T. Kailath ayant pour titre : < ( Detection and estimation in sensor arrays using 
w weighted subspace fitting ) ). 
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ANNEXE A 
GLOSSAIRE 

A.1. Conventions standards : 

x* : conjuguee de x 

x' : transposee du vecteur x 

x 1 " : transposee conjuguee du vecteur x 

X* : transposee de la matrice X 

X 1 : transposee conjuguee de la matrice X 

X* 1 : inverse de X 

X* : pseudo-inverse de X : X # = (x t x)" 1 X t 

tr(X) : trace de X 

E[.] : esperance mathematique 

x, X : valeurs estimees 

x* y : convolution des signaux x et y 

5(t) : distribution de Dirac 
sinc() : fonction sinuscardinal 
® : produit de Kronecker 

A.2. Constantes entieres : 

N :nombre de capteurs composant I'antenne 

M : nombre de trajets 

L : longueur de la reponse impulsionnelle 

K : nombre d'echantillons temporels pour composer I'observation 
spatio-temporelle. 

A.3. Variables scalaires : 

Xn(.) observation sur le n' 6me capteur 
s() : signal emis par la source 
b() : bruit 

a 2 : variance du bruit 

h n () : reponse impulsionnelle du n l6me canal 
Fo : frequence d'echantillonnage 
T e : periode d'echantillonnage 
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am : amplitude du m trajet 

<|> m : dephasage du vn ]hm0 trajet 

9 m : azimut du m i6me trajet 

y m : site du m l6me trajet 

T m : retard de groupe du m l6me trajet 

a n (9) : gain du n i6me capteur dans la direction (9) 

J() : critere a optimiser 

A.4. Variables vectorielles et matricielles : 

x :observation spatiale sur I'antenne, de longueur N 

¥ :observation spatio-temporelle sur I'antenne, de longueur NK 

h : vecteur des reponses impulsionnelles, de longueur N(L+1) 

a(9) : vecteur directionnel spatial, de longueur N 

m(x) : vecteur temporel, de longueur L+1 

l q : matrice identite de taille qxq 

R x ; matrice de covariance du vecteur x 

R y ; matrice de covariance du vecteur y 

R h : matrice de covariance du vecteur h 

M : matrice des M vecteurs temporels : M = [rn(T-|), Q2(tm)] 

A : matrice des M vecteurs directionnels : A = [a(9i), a(9|vi)] 

£2 : matrice des attenuations des M trajets 

H : matrice des filtres mis en colonnes, de taille (L+1)xN 

G : matrice formee des reponses impulsionnelles 
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ANNEXE B 

LE MAXIMUM DE VRAISEMBLANCE DETERMINATE 
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En considerant un signal de la forme : 
y k = A(8)s k + b k 

et le bruit comme etant gaussien, centre de variance a 2 , la 
distribution de probability du bruit est de la forme : 



f(b k ) = 



1 



n 



M 



■exp 



La distribution de probability des donnees qui est a maximiser est 
donnee par la relation : 



1 



f ( x kl s k) = n^TTi ex P 



ML2 

71 a 



-T|Xk-As k | 2 



La vraisemblance des donnees pour le vecteur des donnees 
concatenees est alors : 

K , -J ? | Xk -As k | 2 



L(x 1) x 2l ... t x lc )= n- 



k=17t a 



Pour maximiser la vraisemblance cela est equivalent a minimiser 



la relation : 



-log^x-,, x 2 , .... x k )) = -MKIoga 2 + £|x k - As k | 

k=l 

La minimisation de ce critere par rapport a 6 amene a minimiser : 

K 2 
*K = Zl x k" As k| 

=>s k =(A r A) A f x k 



d'ou : 

K 

e K = Z 
k=1 

comme 



= ZKI 



k=1 



K| 2 = U S = tr ( u k u k) =tr ( u k u k) 



55 
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e k devient : 




Revendications 

1. Precede de localisation spatio-temporelle d'un radio-mobile relativement a une station de base emettrice/r6cep- 
trice caracterise en ce qui! consiste : 

a calculer (1) la r6ponse impulsionnelle du canal de propagation entre le radio-mobile et la station emet- 
trice/r6ceptrice, 

a estimer (2) a partir de la reponse impulsionnelle les couples (angles, retard) du trajet de propagation direct 
et des multitrajets resultats de reflexions sur des obstacles situes dans le champ de propagation des ondes 
entre la station de base et le radio-mobile, 
^ calculer (2) pour chaque trajet les attenuations cij, 

et a calculer (3) la distance d s€parant la station emettrice/r€ceptrice du radio-mobile en consid6rant seule- 
ment les retards et les attenuations du trajet direct et d'un multitrajet quelconque. 

2. Proced6 selon la revendication 1 caracterise en ce qu'il consiste a determiner la distance d par la relation : 



est le rapport des attenuations des propagations sur le trajet direct et le second trajet, 
et \i la partie due a la retlexion de I'onde sur le second trajet. 

3. Proc6de selon la revendication 2 caracterise en ce qu'il consiste pour calculer (2) les attenuations a-j et aq des 
ondes sur le trajet & calculer une matrice diagonale £i def inie par le produit matriciel : 



(t 2 - t i) c 




ou t 2 - ^! est r§cart des temps de propagation entre le trajet direct et le trajet ref lechi, 
c est la cglerite des ondes transmises sur le canal de propagation, 



a 2 



Q=M* HA 1 



ou 




est la matrice des coefficients de la reponse impulsionnelle de canal, 
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H = [h 1l h 2l —.h N ] 

et A * est la matrice transpose des vecteurs directionnels a„ (6 n , y n ) correspondant aux multitrajets. 

4. Proc6d6 seion Tune quelconque des revendications 1 k 3, caract£ris§ en ce qu'il consiste k determiner la r^ponse 
impulsionnelle du canat (1) au moyen d'une antenne multi-capteurs coup 6 a unfiltre causal. k r6ponse impulsion- 
nelle finie. 

5. Proc6d6 selon Tune quelconque des revendications 1 k 4, caract6ris§ en ce qu'il consiste a effectuer une estima- 
tion des angles et des retards en estimant les paramStres temporels k partir de la connarssance des angles 6 m . 

6. Proc£d6 selon Tune quelconque des revendications 1 k 4, caract6ris6 en ce qu'il consiste k effectuer une estima- 
tion des angles et des retards en estimant les param&res spatiaux 8 m k partir des retards x m . 

7. Proc6d£ selon Tune quelconque des revendications 1 k 4, caracteris6 en ce qu'il consiste k effectuer une estima- 
tion conjoints des angles 6 m et des retards t m . 

8. Proc6d§ selon Tune quelconque des revendications 1 a 7, caract£ris§ en ce qu'il consiste k estimer les 
parametres temporels et les directions d'arriv^e des multitrajets par la mise en oeuvre de I'algorithme connu sous 
lenomde << MUSIC )). 

9. Proc6d6 selon Tune quelconque des revendications 1 & 7, caract6ris6 en ce qu'il consiste k estimer les parametres 
temporels et les directions d'arriv£e des multitrajets par la mise en oeuvre de I'algorithme d£terministe connu sous 
la denomination de << maximum de vraisemblance >>. 
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